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Abstract 

By the concurrent use of two different resummation methods, the composite operator formahsm 

and the Dyson-Schwinger equation, we re-examinate the behavior at finite temperature of the 

0(A^)-symmetric A(/9^ model in a generic D-dimensional Euclidean space. In the cases D = 3 

and D = 4, an analysis of the thermal behavior of the renormalized squared mass and coupling 

constant are done for all temperatures. It results that the thermal renormalized squared mass is 

positive and increases monotonically with the temperature. The behavior of the thermal coupling 

constant is quite different in odd or even dimensional space. In D = 3, the thermal coupling 

constant decreases up to a minimum value diferent from zero and then grows up monotonically as 

the temperature increases. In the case D = 4, it is found that the thermal renormalized coupling 

constant tends in the high temperature limit to a constant asymptotic value. Also for general 

D-dimensional Euclidean space, we are able to obtain a formula for the critical temperature of 

the second order phase transition. This formula agrees with previous known values at D = 3 and 

D = A. 
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1 Introduction 



Considerable progress has been done during the last two decades in our understanding of finite 
temperature field theory. A clear account on the foundations of this presently well established 
branch of theoretical physics is done in Refs.|]I[. One of the outstanding questions for the subject 
is the possible existence of a deconfinement phase transition in QCD. It is expected that at suffi- 
ciently high temperatures quarks and gluons are not bounded inside hadrons, but have a behavior 
somehow analogous to the plasma phase, the so called quark-gluon plasma 0. In connection to 
this expected phenomenon it naturally arises questions concerning the nature of the transition 
from the low temperature colorless hadron gas to the high temperature quark-gluon plasma 
In particular, if in fact it is a phase transition, the question of its order is still open. In pure 
Yang-Mills theory it is a point of consensus that there is a phase transition, but if quarks are 
included, the answer seems to be more difficult. Also if the theory involves quarks and gluons, 
another phase transition may occur and its order depends on the number of quark fiavours: the 
chiral symmetry phase transition. The chiral symmetry is spontaneously broken at zero tempera- 
ture and it is restored at some finite temperature. For two massless fiavours it is a second order 
one, while for three fiavours it is a first order phase transition. For more than three fiavours no 
clear conclusion has been obtained about its order, nor even if there is a chiral phase transition. 
Another important point to be noticed in the the case of the chiral symmetry phase transition 
in QCD, is that the order parameter is composite in the fields and perturbative methods fail. 
Non-perturbative methods are mandatories. 
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In any case, "real" QCD is very difficult to deal with and no definite theoretical results for 
many questions have been obtained up to the present moment. So, simpler models still remain 
very useful "laboratories" to try to get some insights about what happens in the real world, 
particularly if long range, non-perturbative effects are involved, as it should be the case for both 
phase transitions mentioned above. For instance, it is well known that the low energy dynamics of 
QCD may be quite well described by the 0(4) model, since there is an isomorphism between the 
0(4) and the SU {2)XSU (2) symmetry groups, the latter being the group of two flavors of massless 
quarks (in the 0{N) model when the vacuum of the model exhibits spontaneous symmetry 
breaking it is known as the linear a model). 

In a previous paper two of the authors derived expressions for the renormalized thermal mass 
and coupling constant in the Xip'^ model in a D-dimensional fiat spacetime at the one-loop approx- 
imation 0. The main results are that the thermal squared mass increases with the temperature, 
while the thermal coupling constant decreases with the temperature. These authors conjectured 
that in D = 3 it could exist a temperature above which the renormalized coupling constant be- 
comes negative. In this regime the system would develop a ffist order phase transition. As we 
will see latter, in the framework of the vector N-component model at a non-perturbative level, the 
answer seems to be negative. In other words, in this context, there is no temperature at which 
the coupling constant changes of sign, and no ffist order phase transition induced by the thermal 
renormalized coupling constant seems to be possible in D = 3. On the other side, in a slightly 
different context, another result that deserves our attention was obtained by Ananos and Svaiter 
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0. These authors studying the {Xip^ + aip^)^ model at finite temperature up the two-loop approx- 
imation, exhibit the existence of the tricritical phenomenon when both, the thermal renormalized 
squared mass and quartic coupling constant become zero for some temperature. 

The purpose of this paper is to get an extension of the above mentioned result ^ to all 
orders of perturbation theory. This may be done considering the O(A^) -symmetric model in the 
framework of the large-N expansion at the leading order in j^, and using resummation methods: 
the composite operator formalism (CJT) and also the Dyson-Schwinger equation, both adapted 
to finite temperature. It is well known that these resummation methods can solve the problem 
of the breakdown of perturbation theory in some massless field theories at finite temperature. 
For instance, in the N-component scalar X^p^ model it is possible, for large N, to sum a class of 
Feynman diagrams, the ring diagrams, by the use of the recurrent Dyson-Schwinger equation. This 
allows to solve at the leading order in j^, the problem of infrared divergences |p. An alternative 
method that takes into account leading and subleading contributions from multiloops diagrams 
was developed by Cornwall, Jackiw and ToumbouUis (CJT formalism) 0. In this approach one 
considers a generalization of the effective action r((/j) which depends not only on the expectation 
value of the field, but also on the expectation value of the time ordered product of two fields. In 
this formalism we naturally sum a large class of diagrams, and the gap equation is easily obtained 
using a variational technique. This formalism has been further developed by Pettini and more 
recently by Amelino-Camelia and Pi 0. It may also be mentioned that recently discussions on 
the pure scalar model have been done by many authors Of special interest are Ref. |11] and 
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Ref.||12||. In this last article Eylal et al. |T2| discussed the phase structure of the 0{N) model in 
a generic D- dimensional euclidean space, and also the sign of the second coefficient of the large-N 



renormalization group (3 function. It may be noticed that Kessler and Neuberger [13| using lattice 



regularization obtained quite surprisingly that for D = 3 the sign of this cofficient is arbitrary. 



The authors in Ref. ||12[ raised the possibility that this ambiguity could be an artifact of the large 
N-expansion. 

In the next sections we apply the above mentioned resummation methods to get non-perturbative 
results for the renormalized thermal squared mass and coupling constant for the vector massive 
X(f'^ model defined in a generic D-dimensional Euclidean space. The outline of the paper is the fol- 
lowing. In section II we briefly discus the CJT formalism. In section III the thermal gap equation 
is obtained. In section IV using the gap equation we sum all the daisy and super-daisy diagrams 
to obtain the thermal renormalized coupling constant. Conclusions are given in section V. In this 
paper we use ^ = c = /c^ = 1. 

2 The Cornwall, Jackiw and Tomboulis (CJT) formalism 



In this section we will briefly discuss the effective action formalism for composite operators 



extended to finite temperature in a D-dimensional Euclidean space |]14[. Let us consider the 
vacuum persistence amplitude Z{J, K) in the presence of local and non-local sources J(x) and 
K{x,y) respectively, where J{x) couples to ^(x) and K{x,y) to |$(x)<I>(y). This object is a 
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generalized generating functional of the Green's functions of the model i.e., 
Z{J,K) = J D$exp|- /($) + J d^xmx)J{x)) + ^J d^x J d^y {<^>{x)K{x,y)<^>{y)) 



, (1) 



where is an appropriate functional measure. Introducing the interaction Lagrange density 
Lint, the classical Euclidean action /($) is given by 



/($) ^Jd'^xJ d^'y {^{x)D^\x - yMy)) + J d'^xLint, 



(2) 



where Dq{x — y) is the free propagator i.e., 



D^\x-y)^-{n + m')5''{x-y). 



(3) 



The generalized effective action r((/7, G) is defined by a double Legendre transform of the gen- 
erahzed generating functional of the connected correlation functions W{lp,G) — \nZ[(p,G) (the 
Helmholtz free energy). 



5W{J, K) 



5J{x 
5W{J, K) 1 



ip{x), 



SK{x) 2 



{ip{x)(p{y) + G{x,y)) 



(4) 
(5) 



and 



r(^, G) = W{J, K)- j d^'x {v{x)J{x)) -\jd''xj d^'y (^(xMy) + G(x, y)) K(x, y). (6) 



In the absence of sources, the generalized effective action r((/?, G) satisfies 

5V{^,G) 



5ip{x) 



J{x)- j d''y{K{x,yMy))^0 



(7) 



and 

5T{^,G) 1 



K{x,y) = ^. (8) 



5G{x,y) 2 

In eq.(§) and eq.(|), the quantity ^{x) is the normahzed vacuum expectation value of the field, 

V9(x) =< 0|$(x)|0 > (9) 

and G{x, y) is the two-point function. If we assume that translation invariance is not spontaneously 
broken we have if{x) = ip and the propagator depends only on the distance in Euclidean space, 
G{x, y) = G{x — y). Consequently it is possible to define the generalized effective potential V{ip, G) 
as a straightforward generalization of the usual definition, 

r(</^, G) = -V{p, G) J d^x. (10) 

The stability conditions given by eqs. (|^) and become in terms of the effective potential, 

dV{v,G) _ dV{v,G) _ 

' dG{k) ' ^ ' 

where G{k) is the Fourier transform of G{x^ y) 

G{x,y) = Gi^-y) = ^j d^ke~^'^--y^Gik) = J ^""^^"^'^^'"'^ (12) 



and = + fyj^. 



2~ 

Let us consider the model described by the Euclidean Lagrange density in the form 



whose effective potential is given by 



^(M^ - - -V2^)G'(a;, x) + ^^(a;, x)G'(x, x) 
2 z o 

Combining the above equation with the stationary requirements given by eq.(|TT|) we have 

^^^ = ^(„2 + ^^^ + ^G(x,.))=0, (15) 

9M2 2 \^ 2^ 2 ^ ' 7 ^ ^ 

The effective potential is obtained by evaluating V{ip,M) from eq. (p!6D . It is composed by the 
classical, the one-loop and two-loop contributions, i.e., 

V{^,M{^)) = Vo + Vi + Vii, (17) 

where the classical contribution is 

Voiip) = \m'^' + A/, (18) 

the one loop contribution is 

Vj{v, M{^)) = \j ^ \n{k' + M\^)\ (19) 
and finally the two-loop contribution is 

V}/((/?, M{ip)) = -^G{x, x)G{x, x). (20) 



In the above expressions we have, 



and G{x, x) is given by 



MH'fi)^m' + ^^' + ^G{x,x), (21) 



To study finite temperature effects there are two different formahsms, the real time formahsm and 
the imaginary time (Matsubara) formahsm. In the Matsubara formalism the Euchdean time r is 
restricted to the interval < r < j3 = ^ and in the functional integral the field $(r, x) satisfies 
periodic boundary conditions in Euclidean time, 

<l>(0,x) = $(/?, x). (23) 

All the Feynman rules are the same as in the zero temperature case, except that the momentum- 
space integral over the zero-th component is replaced by a sum over discrete frequencies. For 
boson fields we have to perform the replacement 

The thermal gap equation may be get from the zero temperature gap equation, 

M=(,).,„^,^,^,^/g,_i_. (25) 



after performing the replacement given by eg . (|2^) . In the next section we use these result to 
analyse the thermal behavior of the squared mass and coupling constant. We remark that the 



previous analysis may also be done using a simpler but ad-hoc procedure |T^, replacing the thermal 
mass obtained in the one-loop approximation by the Dyson-Schwinger equation. Nevertheless, the 
CJT formalism provides a more elegant way and and a consistent basis for this study, since the gap 
equation is derived from a stationary requirement. Indeed, recently Campbell-Smith analysed the 
composite operator formalism in {QED)^ and proved that it reproduces the usual gap equation 
derived using the ad-hoc Dyson-Schwinger approach |jl7 . 



3 The thermal gap equation in the Xcpj^ model 

Let us suppose that our system is in equilibrium with a thermal bath. At the one-loop approxi- 
mation the thermal mass and coupling constant for the Xip"^ model in a D-dimensional Euclidean 
space have been obtained in a previous work and are given by 

= ml + ^ E (^)""V._,(mon/3) (26) 

and 

A(,3) = A„--p-iL^|;(^)" K._,M (27) 

where Ki,{z) is the modified Bessel function and ttiq and Aq are the zero temperature renormalized 
squared mass and coupling constant respectively [^. It is possible to improve the above results 



studying the gap equation for the temperature dependent squared mass. One way to get this 
improvement is to take the 0{N) model, as it has been done by Dolan and Jackiw. In this case 
the Feynman diagrams are classified according to its topology and each of these classes is associated 
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to a given power of In the large limit, only some classes of diagrams, those associated to the 
smallest power of give the leading contributions. These contributions, which include those from 
all daisy and super-daisy diagrams may be summed up, resulting that in the large limit the 
results are exact. This approach was further developed by Weldon and also by Eboli et al. and 



more recently by Drummond et al. |T^. We can obtain from the thermal gap equation resulting 
from the combination of eq.(p4D and eq. (pS]) , an expression to replace the one-loop result above. 
For simplicity we will first analyse the thermal corrections in the the disordered phase, which 
corresponds to absence of spontaneous symmetry breaking. In this case we adopt the notation 
m^{(3) for the squared thermal mass and we suppress the quadratic term in the field. We follow 
a procedure analogous to that Malbouisson and Svaiter have used in p. We use a mix between 
dimensional and zeta function analytical regularizations to evaluate formally the integral over the 
continuous momenta and the summation over the Matsubara frequencies. We get a sum of a polar 
(temperature independent) term plus a thermal analytic correction. The pole is suppressed by 
the renormalization procedure. Then after some technical manipulations, the gap equation may 
be rewritten in the form, 

= ml + E (m(/3)n/3) . (28) 

To solve the above equation, and consequently to go beyond perturbation theory, let us take an 



integral representation of the Bessel function given by 

7r3 1 roo 



2' 
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which is vahd for Re{u) > — ^ and \arg{z)\ < |. This integral representation attend for our 
purposes if we restrict ourselves to D > 1. Substituting eq.(|^) in eq.(|^) and defining 



1 1 



FiD) = ^n-i^o^j^^D^y (30) 
the gap equation becomes 

m'{f3) =ml + XoF{D) {m{P)f-' dt{t' - l)^ ^^J,, _ ^ . (31) 

Defining a new variable r = m{l3)(3t it is easy to show that 

m\i5) = + AoF(D) (m(/?))^-^ / dr[ (— ^)^ - 1 (32) 

When D is odd, the power = p is an integer and the use of the Newton binomial theorem 

will give a very direct way for evaluating m^{(3). When D is even (the most interesting case) the 



D-3 

expansion of ((^^pp)^ l) ^ yields a infinite power series, and the expression for the thermal 



squared mass becomes 



m\(3) = ml + XoP'-^ ^ /(^' k) imiP)/3f' / dr^—, (33) 

^ Jmif3)f3 - 1 



where 



f{D,k) = F{D){-l)'CU, (34) 

2 



and 
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are a generalization of the binomial coeficients. Note that for small values of k the integral that 
appear in eq. (p3|) is a Debye integral of the type 

_,._), (36) 



which is valid for x > and n > 1 [|T^. For k satisfying k > which corresponds to n < 1 in 
the preceeding equation, it is necessary to generalize the Debye integral (the case tt, = is trivial). 



Let us investigate the case n < 0. This generalization has been done by Svaiter and Svaiter |T9 
and the result reads, 

h{x,n)= dr— -— = - J2 -T T-r^Bnlnx + jn^, (37) 

(for odd n), Re{x) > 0,2n > \x\ > and 'jn-i being a constant. The quantites Bn are the Bernoulli 
numbers. Note that this generalization can be used only for high-temperatures i.e. m{l3)(3 < 27r. 
Thus, in the high temperature regime, if we define 



Ii{x,D -3-2k), for x>0, k< 



D-A 



I{x,D -3-2k) = { (38) 

l2{x,D -3-2k), for < x < vr, k>^, 

we may write 

oo 

m\(3) = ml + ^ /(^' ^) i^iP)^f' ^ i^iP)f^, D - 3 - 2k) . (39) 

A:=0 

The above equation gives a non-perturbative expression for the thermal squared mass in the high 
temperature regime in the case of even dimensional Euclidean space. In the odd dimensional case 
the summation in k finishes at 
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For any dimension, it is possible to perform a numerical analysis of the behavior of the renor- 
malized squared mass for all temperatures using eq.(28). It is found that in both cases D = 3, 
and D = 4, the thermal squared mass appears as a positive monotonically increasing function of 
the temperature. For a negative squared mass mp, the model exhibits spontaneous breaking of 
the 0{N) symmetry to 0{N — 1). Since the thermal correction to the squared mass is positive, 
the symmetry is restored at sufficiently high temperature. The critical temperature f3~^ is defined 
as the value of the temperature for which M'^{(3) vanishes. In the neighborhood of the critical 
temperature we can use without loss of generality, eq. (^8D instead of eq . (^51) , since the normalized 
expectation value of the field (the order parameter) vanishes. Using the limiting formula for small 
arguments of the Bessel function it is not difficult to show that the critical temperature in a generic 
D-dimensional Euclidean space {D > 2) is given by: 

(ft- ^ - J^. (40) 

where g{D) = —^T(^ — 1)C(-D — 2). In D = 4 the above result reproduces the known value of 
the critical temperature. For D = 3 the zeta function in g{D) has a pole and a renormalization 
procedure implies that the quantity /^■^^^ is proportional to the the regular part of the analytic 
extension of the zeta function in the neighborhood of the pole. It is not difficult to show that in 
this case the critical temperature is given by = P~^ln{fi^P'^). This result is in agreement with 



the estimate of Einhorn et al. |20|, which has been recently confirmed by numerical simulation 
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4 The thermal couphng constant for the Vector Xipj^ model 

In this section we investigate the behavior of the renormahzed thermal couphng constant of the 
vector N-component Xipjj model. Again, withouth loss of generality, let us suppose that we are 
in the symmetric phase i.e ttIq > 0. To go beyond perturbation theory, we take the leading order 
in in which case we know that the contributions come only from some classes of diagrams 
(the chains of elementary four-point bubbles) and that it is it is possible to perform summations 
over them. Proceeding in that way, we get for the thermal renormahzed coupling constant an 
expression of the form, 

where 

L{m\(iip) = -^(^|:(^)"%_.M/?)n/?), (42) 

and Ao is the zero-temperature renormahzed coupling constant. For simplicity we have suppresed 
the factor j^: everywhere. Incidentally, before going into more details of the non-perturbative finite 
temperature case, some comments are in order. We remark that at zero temperature and D — 4 
it is known that the theory has tachyons, which are related to the occurence of a Landau pole in 
A(m^). The conventional way to circunvect this problem is to interpret the model as an effective 
theory (introducing a ultraviolet cut-off) which restricts the energy to a region far below the 
tachyon mass. This is a necessary procedure, to take into account the widelly spread conjecture 
that the A(/?^ model has a trivial continuum limit at D — 4 and is meaningful only as an effective 
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theory. Coming back to finite temperature, firstly let us investigate the thermal behavior of the 
coupling constant in a Euclidean space satisfying D > 3. Proceeding in a manner analogous as we 
have done in the previous section, we use again an integral representation of the Bessel function 
given by eq.(|2^), which leads to the result, 

L{m' = G{D) {m{P)f-" dt{e - (43) 



where 



2 ' 

Note that there are no poles in the Gamma function and G{D) never vanishes. Defining 



g{D,k) = G{D){-lYC%^ (45) 

2 

it is not dificult to show that 

L (m\(3),p) = (3-^ 9{D, k) {m{(3)pf'^' / dr^—. (46) 
^ ^ Jm(0)0 — 1 



k=0 

A straigthforward calculation gives us the following expression 



L {m'iP),/3) = r E 9iD, k) imiP)/3f'^' I (m(/3)/?, D - 5 - 2k) , (47) 

k=0 

where again I{x, n) is defined in eg. (pH]) . Then, substituting eq.(^7D in the eq.(pD we get the high 
temperature thermal coupling constant for D > 3. In the case D < 3, the integral representation 
of the Bessel function gived by eq.(^) can not be used. Consequently, we take another integral 
representation of the Bessel function i.e. 

K,{z) = -i-] / e~*-tr r^'^+i) (48) 



2 \2J Jo 
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which is vahd for \arg{z)\ < ^ and Re{z^) > 0. A straighforward calculation gives 

L {m\l3),p) = QiD)mif3f-^ dt e'' r^+i (^03(7r,e-^) - l) , (49) 
where the theta function Q3{z,q) is defined by ||18| : 



oo 

,2 



Qs{z, = 1 + 2 ^ cos(2nz). (50) 



n=l 



and Q{D) = ~| (2^)i> • Substituting L {w? {(3) , [3) in eq.(^TD we have a closed expression for the 
thermal renormalized coupling constant in the case D <?>. It is interesting to note that a long time 



ago Braden |22] has employed also the theta function 03(z, 5), to investigate the finite temperature 
\tp^ model. We would like to stress that the behavior of the thermal renormalized coupling 
constant is quite different from the monotonically increasing in temperature behavior obtained for 
the squared mass. Indeed, the kind of the thermal behavior of the coupling constant depends on 
the Euclidean dimension. For D = ?> the coupling constant (as a function of the temperature) 
decrease until some minimum value and then start to increase. This thermal behavior of the 
coupling constant is ploted in fig.(l). For D = 4 the thermal renormalized coupling constant 
tends to a constant value in the high temperature limit. See fig. (2). Our result is consistent with 
the work of Fendley ||23|| . 



It must be noticed that the thermal behaviour of the coupling constant is very sensitive to 
the thermal behaviour of the mass. As an ilustration of this fact we exhibit in fig. (3) the general 
aspect of the coupling constant as function of the temperature for the same model we have treated 



here, but subjected to Wick ordering [Q. In this case all tadpoles are suppressed and the thermal 
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behaviour of the couphng constant does not depend at all on the mass thermal behaviour. We see 
that in this situation the coupling constant is a monotonic decreasing function of the temperature. 
The absence of Wick ordering deeply changes this behaviour. We show in fig. (4) for D = 3 in 
the same scale the plots for X{T) with and without Wick ordering, respectivelly the lower and 
the upper curves. In the region of temperatures where Xw{T) goes practically to zero, A(T) is 
practically constant at a value slightly lower than the common zero-temperature coupling constant 
Aq. The growth of A(r) with the temperature presented in fig.(l) is in a much smaller scale for 
A(T) than in fig. (4). In fact this growth is "microscopic" in a scale where the Wick ordered 
coupling constant Xw{T) presents asymptotic thermal freedom. 

5 Conclusions 

We have done in this paper an analysis of the vector X(p'^ model in a flat D-dimensional Euclidean 
space in equilibrium with a thermal bath. The form of the thermal corrections to the mass 
and coupling constant have been discussed using resummation methods. We have formulated 
the problem in a general framework, but our results are at leading order in the expansion. 
We have chosen this way of working, in order to get answers as much as possible of a non- 
perturbative character. In what concerns the thermal mass behavior, we have shown that the 
thermal renormalized squared mass is a monotonic increasing function of the temperature for any 
Euclidean dimension. We have been able to obtain a general formula for the critical temperature 
of the second order phase transition valid for any Euclidean dimension D > 2, provided the 
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necessary renormalization procedure is done to circunvect singularities of the zeta function. The 
values obtained for the critical temperatures in D = 3 and D = 4 agree with previous results. 

The behavior of the thermal coupling constant depends on the dimensionality of the Euclidean 
space. In D = 3 the renormalized coupling constant decreases until some positive minimum 
value and then starts to increase slightly as a function of the temperature. See fig.(l). This 
result seems to indicate that at a non-perturbative level (for large N) in the framework of the 
vector N-component model, the answer to the question raised in Ref.(|^) is negative: there is no 
first order phase transition induced by the thermal coupling constant in D = 3. In D = 4 the 
thermal renormalized coupling constant in the high temperature limit tends to a constant value, 

r- 3 



(Ao — (3-^)Ao, which coincides exactly with the result obtained by Fendley ||23[| ). 



A natural extension of this work should be to go beyond the leading order results using 
renormalization group methods. Another possible direction is to introduce an abelian gauge field 
coupled to the N-component scalar field. In D = 3 a topological Chern-Simons term may be 
added, and also a tp^ term. In this case we have shown in a previous work using perturbative 
and semi-classical techniques that the topological mass makes appear a richer phase structure 
introducing the possibility of first or second order phase transitions depending on the value of the 



topological mass The use of resummation methods to investigate the thermal behavior of 

the physical quantities could generalize these previous results to the N-component vector model. 
These will be subjects of future investigation. 
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Figure 1: Coupling constant thermal behavior obtained from eqs. (p8D ,(^lD and (^21) in dimension 
D = 3 . 
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Figure 2: Coupling constant thermal behavior obtained from eqs. (|28D ,(HlD and (^21) in dimension 
D = A. 
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Figure 4: Compared thermal behaviors of the couphng constant for the Wick-ordered model and 
non Wick-ordered model. 



